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Experimental tests of QCD through its predictions for the strange-quark content of the proton 
have been drastically restricted by our lack of knowledge of the violation of charge symmetry (CSV). 

We find unexpectedly tiny CSV in the proton’s electromagnetic form factors by performing the first 
extraction of these quantities based on an analysis of lattice QCD data. The resulting values are an 
order of magnitude smaller than current bounds on proton strangeness from parity violating electron- 
proton scattering experiments. This result paves the way for a new generation of experimental 
measurements of the proton’s strange form factors to challenge the predictions of QCD. 

PACS numbers: 13.40.Gp, 12.39.Fe, 14.20.Dh 


Charge symmetry is the invariance of the strong in¬ 
teraction under an isospin rotation exchanging u and d 
quarks (i.e., exchanging the proton and neutron). The 
violation of this symmetry (CSV) is arguably small: the 
proton-neutron mass difference is one part in a thou¬ 
sand [1] and many nuclear reactions proceed identically 
if protons and neutrons are interchanged. The effects of 
this small CSV, however, may be hugely significant. For 
example, if the proton-neutron mass difference were re¬ 
versed protons could decay and atoms could not form. 
Charge symmetry violation also explains the discrep¬ 
ancy between the calculated and measured binding en¬ 
ergy differences of mirror nuclei (Okamoto-Nolen-Schiffer 
anomaly) [H [5] and may play a role in precision tests of 
the Standard Model [1], including those at the LHC [5]. 

In the late 1980s it was suggested that one could use 
measurements of neutral weak current matrix elements 
by parity-violating electron scattering (PVES) [MS] to 
determine the contribution of strange quark-antiquark 
pairs to the elastic electroweak form factors of the nu¬ 
cleon. These ‘strange form factors’ have been the focus of 
intensive experimental and theoretical effort for the past 
two decades [5]. At present, the accuracy of theoretical 
calculations of these quantities [TT)UTS] exceeds that of 
the best experimental values [14] by almost an order of 
magnitude—a remarkable exception in strong-interaction 
physics. The limiting factor in future state-of-the-art 
PVES measurements at Mainz [m ns] and JLab [nh 
m is theoretical, arising from the assumption that CSV 
in the proton’s electromagnetic form factors is negligible. 

Precisely, CSV form factors Gcsvj if not accounted 


for, mimic the strange-quark contribution in the 

combination of form factors accessed by experiment: the 
measured neutral weak current matrix elements G^^ 
may be expressed as 

= (1 - 4sin^ dw) - G^’/V - G|;/m + Gcsv, 

( 1 ) 

where the weak mixing-angle, 9w, and the total elec¬ 
tromagnetic form factors, , are precisely deter¬ 

mined from other experimental studies. With theoret¬ 
ical predictions of the size of Gcsv varying through sev¬ 
eral orders of magnitude [201122] . this uncertainty has 
halted experimental parity-violating electron scattering 
programs [TU] . 

In this Letter we report the hrst determination of CSV 
in the proton’s electromagnetic form factors based on an 
analysis of lattice QCD data. In terms of individual u and 
d-quark contributions to the Sachs electric and magnetic 
form factors of the proton and neutron (conventionally 
defined without the charge factors), the CSV form factors 
which we calculate are dehned as 

XU _ ^P,u ^n,d rd _ r^P,d ^n,u /,y\ 

°E/M — ^E/M ^E/M^ °E/M “ ^EjM ^E/M^ G) 

where we explicitly calculate and perform the 

subtractions indicated. The combination relevant to ex¬ 
perimental determinations of nucleon strangeness using 
Eq. Q is 
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The lattice results used here are an extension of those 
reported in Refs. [23l [24]; we include two independent 
sets of 2 + 1-flavor simulations at different values of the fi¬ 
nite lattice spacing a. Any discretization artifacts should 
appear at 0{a^). Each set consists of results for the indi¬ 
vidual connected quark contributions to the electromag¬ 
netic form factors of the entire outer-ring baryon octet at 
a range of pion masses down to 220 MeV and at 6 (set I) 
or 7 (set II) fixed values of the momentum transfer up 
to 1.4 GeV^. These values of are relevant to experi¬ 
mental studies of the strange nucleon form factors. The 
lattice volumes are L^xT = 32^ x64 and 48^ x96, and the 
lattice spacings are a? = 0.0055(3) fm^ and 0.0038(2) fm^ 
(set using various singlet quantities [25l |26|) for the two 
sets respectively. 

Our extraction of the CSV form factors from the lat¬ 
tice simulations is based on the extrapolation of those re¬ 
sults to infinite volume and to the physical pseudoscalar 
masses using a formalism based on connected chiral per¬ 
turbation theory 113 HHj- The extrapolation procedure 
is detailed in Refs. ^51 [H] . The small finite-volume cor¬ 
rections are model-independent and the chiral extrapola¬ 
tion is demonstrated to be under control—the fit includes 
lattice data at low meson masses within the convergence 
regime of the effective theory, and it reproduces the ex¬ 
perimental form factors at the physical masses [H] |13| . 
To determine the CSV terms we must extend that work 
to incorporate the breaking of the flavor-SU(2) symme¬ 


try, i.e., to allow for unequal light quark masses, m„ ^ 
md- This is a simple extension, and is performed precisely 
as in previous work where the same procedure was used 
to evaluate the mass splittings among members of baryon 
isospin multiplets jUj , the CSV sigma terms HO] , and the 
CSV parton distribution moments [3T| from 2+1-flavor 
lattice simulation results. In brief, the low-energy param¬ 
eters which appear in the SU(2)-breaking terms in the 
chiral extrapolation expressions for the CSV form factors 
also appear in the isospin-averaged expressions. These 
parameters are thus fixed by the hts to the Nf = 2 + 1 
lattice QCD simulations on the baryon octet which are 
presented in Refs. [2^ 124] . 

In principle, the CSV form factors on an infinite vol¬ 
ume and at the physical pseudoscalar masses may thus, 
given the extrapolations of Refs. [531 lU) be evaluated 
simply by performing the subtractions shown in Eq. Q. 
This procedure, however, suffers from a significant sys¬ 
tematic effect resulting from the omission of quark-line 
disconnected contributions in the simulations. To ac¬ 
count for this omission we use the chiral extrapolation 
expressions to model the disconnected pieces of the loop 
integral expressions. This amounts to the replacement of 
the ‘connected’ extrapolation coefficients of Refs. [231124] 
by the ‘full’ expressions, where the free parameters re¬ 
main as fixed by the connected hts. The resulting ex¬ 
pressions for the CSV electric and magnetic form fac¬ 
tors (including disconnected quark-line contributions) as 
a function of meson masses can be written as 
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where all symbols, including the low-energy constants cf^", are dehned in Refs. [531 HI]- The leading-order 
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FIG. 1. Diagrammatic quark-line skeleton representation of 
omitted contributions to the CSV form factors. Solid and 
wavy lines represent quarks and photons respectively. The 
crosses denote quark mass insertions, i.e., the figures repre¬ 
sent the contribution from disconnected quark loops to CSV 
arising from the different {u and d quark) masses of: (a): the 
struck sea quark; (b): spectator quarks. These contributions 
are proportional to BijUd — mu)- 


loop integral expressions include meson loops with octet- 
baryon [Iq) or decuplet-baryon {Id) intermediate states, 
as well as tadpole loops {It)- The Gell-Mann-Oakes- 
Renner relation suggests the definition 

B{md - mu) = ( 8 ) 

where R denotes the light-quark mass ratio R = mu/md- 
We take R = 0.553(43), determined by a fit to meson 
decay rates (32] . The final results are all consistent within 
uncertainties if we instead take the FLAG value R = 
0.46(2)(2) [33]. 

All of the low-energy parameters, other than 

and , are determined from the chiral hts 

to the connected contribution to the isospin-averaged 
electromagnetic form factors which are described in 
Refs. [331121] ■ While this procedure systematically in¬ 
cludes some of the disconnected contribution to the GSV 
form factors, other disconnected terms—those which are 
linear in B{md — mu) and not generated by chiral loga¬ 
rithms from meson loops—cannot be determined in this 
way. Precisely, the terms which are generated by the La- 

.^1 rc- ■ ^ E/M E/M J E/M 

grangian pieces with coemcients and 

cannot be determined from the present lattice simula¬ 
tions. Physically, these terms arise from the diagrams 
illustrated and described in Fig. These contributions 
are anticipated to be small based on the success of va¬ 
lence quark models in reproducing form factor data. This 
is also supported by the results of direct lattice QGD 
calculations of Ge/m which find that the disconnected 
contributions at small finite momentum transfer are con¬ 
sistent with zero and are bounded at the 1% level [51] . 
The terms corresponding to the low-energy parameters 
and are only part of that small dis¬ 

connected contribution. 

We choose to set contributions from the unknown 
and terms to 0, with an uncertainty 




FIG. 2. Quark-line skeleton diagrams of the meson loops used 
to model the omitted contributions to the CSV form factors. 
Solid and wavy lines represent quarks and photons respec¬ 
tively. The crosses denote quark mass insertions into: (a): 
the struck sea quark in the meson loop; (b): the meson loop 
spectator quark. 


taken to be twice the magnitude of the corresponding 
contributions from meson loop diagrams, evaluated with 
a dipole cutoff regulator with mass scale A = 0.8(2) GeV. 
We suggest that this error estimate is extremely conser¬ 
vative. The use of this method to evaluate the loops is 
justified by the well-established and successful use of this 
model to relate full and partially-quenched lattice QGD 
calculations [35]. The loop diagra m used to estimate the 
0 ^ 2 ^ terms is represented in Fig. 2(b) where only the 
‘loop spectator’ quark mass (i.e., the valence-quark part 
of the meson mass) is changed. For the term, rep¬ 

resented in Fig. |2(a)[ only the sea-quark part of the loop 
meson mass is considered. These contributions are added 
in quadrature. The magnitude of this contribution to the 
total uncertainty varies with it is largest at our low¬ 
est Q^-values where it contributes 20-60% of the quoted 
uncertainty on the final result (depending which of 
one is considering), while at larger values of consis¬ 
tent with the suppression of meson loops at high-Q^, it 
contributes 1-15%. 


The results of this analysis for the individual u and 
d-quark contributions to the CSV electric and magnetic 
form factors of the proton are shown in Fig. The close 
agreement of the two sets of simulations (at different 
lattice spacings a and on different simulation volumes) 
confirms that the finite-volume corrections and chiral ex¬ 
trapolations are under control and that any discretization 
effects resulting from the hnite lattice spacing are small. 
The size of the CSV form factor combination, GcsVi rel¬ 
evant to PVES experiments probing the strange electric 
and magnetic form factors of the nucleon by Eq. Q, 
is shown in Fig. [^ This result gives quantitative con¬ 
firmation that CSV effects in the electromagnetic form 
factors, for momentum transfers up to approximately 
1.4 GeV^, are at the level of 0.2% of the relevant pro¬ 
ton form factors—an order of magnitude smaller than 
the precision of existing PVES studies. To put this in 
perspective, the level of CSV shown in Fig. |4(b)| is equiv¬ 
alent to a CSV difference in charge radii of less than one 
attometer. These precise results open the door for a new 
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generation of experiments to probe the structure of the 
quantum vacuum through the strange quark form factors. 
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FIG. 3. Individual up and down quark contributions to the CSV form factors. These terms are combined to give the total CSV 
form factors Gcsv = Blue points and green crosses show the results of data sets I and II extrapolated to 

the physical point, with corrections applied to model the omitted disconnected terms. 
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FIG. 4. Magnetic and electric CSV form factors as relevant to experimental determinations of nucleon strangeness. The blue 
circles and green crosses denote our results based on simulation sets I (a^ = 0.0055(3) fm^) and II (a^ = 0.0038(2) fm^), 
respectively. 




































